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Abstract. - We theoretically analyze non-local electron transport in multi-terminal normal- metal- 
superconductor-normal-metal (NSN) devices in the presence of an external ac voltage bias. Our 
analysis reveals a number of interesting effects, such as, e.g., photon- assisted violation of balance 
between crossed Andreev reflection (CAR) and elastic cotunneling (EC). We demonstrate that at 
sufficiently small (typically subgap) frequencies of an external ac signal and at low temperatures 
the non-local conductance of the NSN device turns negative implying that in this regime CAR 
contribution to the non-local current dominates over that of EC. Our predictions can be directly 
tested in future experiments. 



Introduction. — The phenomenon of non-local (or 
crossed) Andreev reflection [1, 2] is known to occur 
in multi-terminal hybrid normal-metal-superconductor- 
normal-metal (NSN) proximity structures and involves 
two subgap electrons entering a superconductor from two 
different normal terminals and forming a Cooper pair 
there. Such crossed Andreev reflection (CAR) manifests 
itself, e.g., in the dependence of the current I\ through 
the left NS interface of an NSN structure on the volt- 
age Vi across the right NS interface. As a result, the 
non-local conductance G12 = —dI\jdVi of an NSN device 
differs from zero and can be detected experimentally. Var- 
ious aspects of this intriguing phenomenon have recently 
become a subject of intensive investigations both in ex- 
periment [3-6] and in theory, see, e.g., [7-19] and further 
references therein. 

It is important that CAR is not the only process which 
contributes to the non-local conductance Gn- Another 
relevant process is direct electron transfer between two 
normal terminals through the superconductor. In the tun- 
neling limit this process is just the well known elastic co- 
tunneling (EC). It was demonstrated [7] that in the lowest 
order in tunneling the contributions from EC and CAR to 
G12 exactly cancel each other in the limit of low temper- 
atures and voltages, i.e. the non-local conductance G12 
tends to zero in this limit. In certain cases this observa- 
tion might complicate experimental identification of CAR 
in NSN structures with weakly transmitting interfaces. 

Despite such possible complications in experiments [3-6] 



the non-local signal G12 was successfully detected both at 
high and low temperatures showing a rich structure of 
non-trivial features. Some of these features are currently 
not yet fully understood and are still waiting for their 
adequate theoretical interpretation. These observations - 
along with various theoretical predictions - also demon- 
strate that the exact cancellation between EC and CAR 
contributions [7] can be violated in a number of ways. One 
of them is simply to lift the spin degeneracy in the prob- 
lem, e.g., by considering NSN structures with spin-active 
interfaces [14] or by using ferromagnets (F) instead of nor- 
mal electrodes [9] . Experiments with FSF structures [3] di- 
rectly demonstrated the dependence of the non-local con- 
ductance G12 on the polarization of F-electrodes. 

In a spin degenerate case the non-local conductance 
G12 does not vanish beyond the tunneling limit [12,14], 
i.e. the cancellation between EC and CAR terms is effec- 
tively eliminated due to higher order electron tunneling 
processes which become significant at higher barrier trans- 
missions. This theory predicts positive non-local conduc- 
tance implying that direct electron transfer should always 
dominate over CAR at higher barrier transmissions. Fur- 
thermore, for ballistic NSN structures CAR was predicted 
to vanish completely in the limit of fully open NS inter- 
faces [12,14]. 

It is worth pointing out that both positive and nega- 
tive non-local signals have been detected in multi-terminal 
NSN devices [3,4,6]. It was argued [15] that negative non- 
local currents could possibly be attributed to the effect of 
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electron-electron interactions. In the presence of strong 
Coulomb interaction negative non-local conductance was 
obtained in single-level quantum dots coupled to normal 
and superconducting electrodes [18]. In general, this im- 
portant issue deserves a detailed theoretical investigation 
which should account for non-trivial interplay between dis- 
order and Coulomb interaction. Previously a similar anal- 
ysis was developed for (local) Andreev reflection in NS hy- 
brid structures [20-22]. This analysis revealed a number 
of interesting features which can also significantly affect 
non-local properties of NSN devices. 

Leaving this important issue for future investigations, 
here we address a somewhat different situation. Namely, 
we will study both EC and CAR processes in NSN struc- 
tures in the presence of an external ac bias. We will 
demonstrate that application of an ac electromagnetic field 
to such structures lifts exact cancellation between EC and 
CAR contributions already in the lowest order in barrier 
transmissions. Under certain conditions CAR can domi- 
nate over EC, in which case the non-local conductance of 
the system G\2 turns negative. 

Non-local currents in NSN devices. — We will an- 
alyze the behavior of multi-terminal NSN structures. An 
example of such a structure is schematically shown in Fig. 
1: A long superconducting wire with BCS order parame- 
ter A(T) (connected to bulk superconducting terminals) is 
coupled to two normal metals via tunnel barriers with re- 
sistances Ri and i?2 ■ We also assume that time dependent 
voltages V\ and V 2 are applied to normal metal terminals. 

The Hamiltonian of this system can be expressed in the 
form 



H = Hi + H2 + Hs + Hi t + H 2 r, 



(1) 



where 



H r = ( 6k ~ eV r( t )) C l,k,a C r,k,a> ( 2 ) 

a=T,l k 

are the Hamiltonians of the normal leads (r = 1, 2), 
k 

(3) 

is the Hamiltonian of the superconducting electrode and 

H r ,T = H [t^U^P + CC. ] (4) 

are the tunnel Hamiltonians describing electron transfer 
between normal terminals and a superconductor. In what 
follows we will restrict our analysis to the spin-degenerate 
case and will also ignore inelastic effects. 

In order to evaluate the non-local current response to 
the applied voltages we proceed within the standard per- 
turbative approach in the tunnel Hamiltonians Hi t and 
i?2,T combined with the Keldysh technique. The whole 
calculation is analogous to that performed in Ref. [7] with 
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Fig. 1: Schematics of the system under consideration. A long 
superconducting wire is coupled to two normal metal leads via 
the two tunnel junctions with the resistances Ri and ife- 



the only important difference that here the applied volt- 
ages Vi,2 explicitly depend on time. Calculating the cur- 
rent I\ across the first barrier perturbatively in the corre- 
sponding tunneling amplitudes, we obtain 



h{V u V 2 ) = /j 0) (Ui) +I n (Vi) + h 2 (V 2 ), 



where the first two terms lf\Vi) and hi(Vi) represent 
the standard contributions to the current across NS inter- 
face evaluated respectively in the first and second orders 
in the interface transmission, while the term I\ 2 {V 2 ) de- 
scribes the non-local contribution to the current across the 
first junction due to the presence of the voltage bias V 2 at 
the second normal terminal. For these three contributions 
we find 



(5) 



(o) _ 



n 
in 



I tr [ - GihGsAA + Gg* t 1 GitiA] t (6) 



tr [ - GxhGsiXGxhGsilk 
■GsAGxhGsAGxiiK] tt , 



I12 = 2 tr [ ~ GihGshGAGsiik 



G s i 2 G2t 2 GsV 1 Gihk] 



(7) 



(8) 



Eqs. (5)-(8) express the current in terms of unperturbed 
Green-Keldysh functions of normal and superconducting 
terminals, respectively G1.2, Gs- As usually, these func- 
tions are 4x4 matrices in Keldysh and Nambu spaces. In 
equilibrium they depend only on the time difference. For 
instance, for the Green-Keldysh function of the supercon- 
ducting terminal one has 

/dE 
— e- lE ^- t ^G s {E,x l ,X2), 



where 



G s (E,x 1 ,x 2 )=G%(E,x 1 ,x 2 ) 



i + Q(E) 



■Gi{E, Xl ,x 2 )® 1 ^ E h z , (9) 



er, is the Pauli matrix and the 2x2 matrices 
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Gg' A (E, xi, X2) read 



fit*) 



Fi 



E 



A iXn(jl)4( a! ') 

(BiiO^-A 2 -^ 



A 3Xn(xi)Xn(x 2 ) 



-£ 2 



(g-4„±iO)x„(x 1 ) x *(x 2 ) 



Here we have denned the normal state electron wave func- 
tions in the superconducting electrode Xn {x) , and the ma- 
trix 



Q(E) = 



1 - 2n(E) 2n(E) 

2 - 2n(E) -1 + 2n{E) 



(10) 



where n(E) = 1/(1 + exp(E/T)) is the Fermi function. 
The matrix Q(E) satisfies the normalization condition 
Q 2 {E) = 1. The Green-Keldysh functions for N-terminals 
G lj2 are defined analogously, one should just set A = 0. 

Let us also note that the trace operation in Eqs. (6)-(8) 
includes taking the trace in 4 x 4 Kcldysh-Nambu space to- 
gether with the convolution over intermediate coordinates 
and times. The subscript t ,t implies that the outer times 
in the products of the Green functions should be equal 
and no integration over these times should be performed. 
Eqs. (6)-(8) also involve the 4x4 matrices A and The 
matrix A is diagonal with the following matrix elements: 
An = A44 = —1, A22 = A33 = 1. The matrices t\i read 



where 









tje 



(11) 



Pl,2(*) 



J dt'eV ia {t') 



are the time dependent phases across the tunnel barriers. 

The second order contributions to the current (7), (8) 
describe several different physical effects. The differ- 
ence between these effects can be illustrated, e.g., in the 
course of averaging of the Green functions over disorder. 
For instance, the term (7) contains the product of four 
Green functions G\GsG\Gs- While averaging, one can 
perform the perturbation theory in 1/gs, l/ffJV; where 
gs,N — 2ir/e 2 rs,N are dimensionless conductances of su- 
perconducting and normal leads (see Fig. 1). Keeping the 
leading and next to the leading contributions, one finds 



(G 1 G S G 1 G S ) = A Q + A 1 + A 2 , 



(12) 



where 



A = (G 1 )(Gs)(G 1 )(G s ) (13) 
is the leading contribution oc (l/gs)°, (1/gjv) , and 



A 1 cx {G 1 G 1 )(Gs)(G s )ocl/g N , 
A 2 cx (G 1 )(G 1 )(G s G s )^l/gs 



(14) 



are the first order corrections. Here the (GjGj) = 
(GjGj) — (Gj)(Gj) are irreducible averages which, in gen- 
eral, contain both diffusons and Cooperons. The aver- 
age A yields the standard subgap current [23] which can 
be disregarded in the limit of small barrier transmissions 
considered here. The contribution A\ describes Andreev 
reflection enhanced by disorder in the N-metal [20,24,25]. 
This contribution is also omitted here. For the local cur- 
rent In this approximation is justified (1) at energies (e.g. 
T or u>) well above the Thouless energy of the N-terminal 
or (ii) for sufficiently strong dephasing in the N-metal 
or (Hi) for g N » g s . Furthermore, this approximation 
docs not affect the lowest order in tunneling contribution 
cx I/R1R2 to the non-local current I12 at all. Thus, below 
we will keep only the terms cx A2 in Eq. (7) and retain 

averages of the type (Gi)(GsGs)(G2) in Eq. (8). 

Evaluating the traces in Eqs. (6)- (8) and assuming that 
exact particle-hole symmetry remains preserved, we get 

= --± f -jdl>G?{t-l!,xux{)h{t-1>) 



eN s R\ 
x sin[v?i(i) - V\{t% 
i 



(15) 



dtidt 2 



2e 3 N 2 R 2 

x [G§(t - fi,xi,xi)Gs(i 2 - t, xi,xi) 
- Fg(t - 1 1 , xi , xi )F A (t 2 - t , xi , xx )] 
x n(ti - t 2 ) sin [ipxih) - ^(fe)] , (16) 

dt\dt2 



2e 3 N 2 R 1 R 2 
x [G§(t - h,x 1 ,x 2 )G A (t 2 - t, x 2 ,xi) 
- Fg(t - h,x u x 2 )Fg(t2 - t, x 2 , xi)] 
x n(ii - t 2 ) sin [<p 2 (h) - ^2^2)], (17) 
where we defined 

*(*)=/ d ^e- iE \l-2n(E)). 

Next we assume that all the electrodes in our system 
are diffusive. In this case we need to average the above 
expressions over disorder. This averaging is accomplished 
with the aid of the following rules 

(Gf(E,x 1 ,x 1 )) = -mN s 

(18) 



VE 2 - A 2 



VA 2 - E 2 



{G§(E u x 1 ,x 2 )G£(E 2 ,X2,x 1 )) = 

S \ W(E 1 )W*(E 2 ), 
xV s (W(E 1 )-W*(E 2 ),x 1 ,x 2 ), (19) 

(Fg(E 1 ,x 1 ,x 2 )Fg(E2,X2,x 1 ))=irN s 
A 2 



W(E 1 )W*(E 2 ) 



V S (W(E 1 )- W*{E 2 ), Xl ,X2pO) 
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Here N$ is the normal density of states in the supercon- 
ducting electrode, 



W(E) = signE9(\E\ - A)y/ E 2 - A 2 

+ i6(A- \E\)VA 2 -E 2 , (21) 

and T>s(u), x.x') is the diffuson, which satisfies the equa- 
tion 



(-iuu - D S X7 2 ) V s (u>, x, x) = S(x - x'), 



(22) 



where Ds is the diffusion coefficient in the supercon- 
ducting electrode. Here we assume that the time rever- 
sal symmetry is maintained in our problem, therefore we 
do not need to distinguish between the diffuson and the 
Cooperon. Then we obtain 



7(0) . . 



/ dt' (G?(t-t',x,x))h(t-t') 
sRi J 



eN s R 
x sin[pi(t) - ipi{t')], 



(23) 



TTl 



2e 3 N s R 2 



dtidtz 



dE 1 dE 2 dE 3 
(2wf 



x e 



hi = 



x e 



iE 1 (t-t- L )-iE 2 (t 2 -i)-iE 3 {ti-t 2 ) 
. EiE 2 - A 2 \ 

xV s (W(E 1 )-W*(E 2 ),x 1 ,x 1 ) 
x em [p^h) (24) 

m f ^ ^ f dE x dE 2 dE z 

2eSN sRl R 2 J dhdh J (27r)3 

iE 1 (t-t 1 )-iE 2 (t2-t)-iE 3 (t 1 -t 2 ) 

E X E 2 -A 2 \ 

x 1 1 + wmw*{E 2) ) n ^ 

xV s (W(E 1 )-W*(E 2 ),x 1 ,x 2 ) 

x sin [ip 2 (ti) - <p 2 (t 2 )]. (25) 

Together with Eq. (5) these general expressions define 
the net current across the first barrier for arbitrary de- 
pendence of the applied voltages V\ i2 (t) on time. Below 
we will analyze these expressions in several specific limits. 

Linear ac response. — Let us choose the applied 
voltages in the form Vj(t) = Vj + 5Vje~ lLOjt . Accordingly, 
the time dependent phases are defined as 



ipj(t) = eVjt + eSV 3 



1 _ e -«^t 



ILUj 



Substituting these expressions into the above results for 
the current and expanding in the amplitudes of the ac 
signal SVj, we arrive at the correction to the current due 
to the presence of ac bias: 

Sh = Gn(wi, Vi^Vie-*" 1 * - Gi 2 (w2, V 2 )5V 2 e-^ 2t . (26) 

Here the choice of the minus sign in front of the second 
term is just a matter of convention which ensures positive 
values of G\ 2 at high bias voltages. 



Let us first specify the nonlocal conductance G\ 2 (oj, V 2 ). 
ft reads 



Gi 2 = / dE I 1 



E 2 - <£L - A 2 



W{E+^)W* (E-%) J 

Vs (W(E+f>-W* (E-%), Xl ,x 2 ) 
X 8e 2 N s RiR 2 LJ 

x n (E - | + eV 2 ) + n (e - | - eV 2 ^j 

-n(E+^ + eV 2 )-n(E+^-r] 2 

The local ac conductance is given by the sum of two terms 
G u (u, V) = G ( V(u, V) + 5Gu(u, V), (28) 

where 

dE 



(27) 



G " = / 



2R l u 



\E\6(\E\ - A) _.E6(A-\E\) 



VE 2 - A 2 VA 2 - E 2 

x [n(E -oj- eVi) + n(E - u + eVi) 

-n(E-eVi)-n(E + eVi)], (29) 



and the correction SGh(uj, V\) is given by Eq. (27) with 
a simple replacement T>$ (w, X\, x 2 ) — » T>g (u>, X\, X\). 

Eqs. (27)- (28) apply for any sample geometry. Let us 
now specify these equations for the NSN structure illus- 
trated in Fig. 1. Provided the superconductor can be 
treated as a sufficiently thin wire with cross-section s, the 
diffuson takes the form 



T> s (uj,x,x') 



cosh 




\x-x'\) 


2sD s 


sinh 


V u s 



cosh 



^f(L-x-x') 



2sD s J^fsmh 



I —ZUJ T 



(30) 



Substituting this expression into Eqs. (27)-(28) we arrive 
at our final results. 

In the zero frequency limit lo — > the above general 
expressions reduce to 



\E\6{\E\ - A) 



dE 



VE 2 - A 2 ^cosh 2 ^ 1 



2 - tanh - tanh *±^'l , (31) 



2r L R\ 
Gi 2 (0,y 2 ) 



2T 



2r L R x R 2 



2 — tanh 
- tanh 



2T 



eV 2 



2T 
A + eV 2 
2T 



(32) 



where the resistances Tl, t Xi , Tl-xi and ri_ l2 are the 
normal state resistances of the corresponding segments of 
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Fig. 2: Real and imaginary parts of the non-local conductance 
Gi2(w,V2) (respectively left and right panels) at zero temper- 
ature and at uj = 0.5A. Here we fix L — 10£, xi = 3£, 
X2 = 4£, where £ = ^/Ds/A is the superconducting coher- 
ence length. The the conductance values are normalized by 
Go = r xl r L - X2 /r L R 1 R 2 . 



the superconducting wire as it is illustrated in Fig. 1. 
We observe that at T — > and eVi < A the non-local 
conductance G12 (32) tends to zero in agreement with [7]. 

At non-zero frequencies the non-local conductance 
G\2(u>,V2) was evaluated numerically. The corresponding 
results are presented in Figs. 2 and 3. Fig. 2 illustrates 
the behavior of the non-local conductance for a given value 
of oj = 0.5A, T — > and different voltage bias values V2. 
We observe that an external ac bias lifts the balance be- 
tween EC and CAR processes essentially at all values V2. 
On top of that, at not very large bias voltages V2 the real 
part of the non-local conductance G12 becomes negative 
implying that in this regime crossed Andreev reflection 
dominates over the elastic cotunneling. At eV<i = ±A =p oj 
this effect reaches its maximum and G12 shows clear dips 
related to the divergence of the superconducting density of 
states at the gap energy. For |eV2 1 > A — oj the real part of 
the non-local conductance starts growing with increasing 
bias and eventually becomes positive at \e\2\ > A — u>/2. 
At eV2 — A the non-local conductance ReGi2 has a peak 
which is again related to the behavior of superconducting 
density of states in the vicinity of the gap energy A. Also 
the behavior of the imaginary part of the non-local con- 
ductance ImGi2 (which differs from zero at |eV^| > A — oj) 
turns out to be rather non-trivial, see Fig. 2. 

The non-local conductance G12 at zero dc bias V2 = 
and different frequencies oj of an ac field is displayed in 
Fig. 3. We again observe that at T = and at frequencies 
oj below the superconducting gap A CAR dominates over 
EC thus turning the real part part of the non-local con- 
ductance negative, see Fig. 3a. This picture gets modified 
at non-zero T, in which case EC contribution turns out 
to exceed that of CAR at sufficiently small oj, though at 
higher frequencies CAR can still win over EC, cf. Fig. 3c. 
We also observe that a peak in the real part of the non- 
local conductance develops at low frequencies and non- 
zero T. The height of this peak is essentially determined 
by elastic cotunneling and is defined by a simple combina- 
tion G [l - tanh A/2T] (sec Eq. (32) j, while its width is 
roughly proportional to the Thouless energy Ds/\x\ — X2I 2 
and slightly increases with temperature. At high enough 




Fig. 3: The same as in Fig. 2 at zero bias voltage V2 = and 
different ac bias frequencies oj: (a) Re Gi2(w,0) at T — 0, (b) 
Im Gi 2 (w,0) at T = 0, (c) Re Gi 2 (w,0) at T = 0.3A , where 
Ao is the value of the superconducting gap at T = 0, (d) Im 
Gi 2 (w, 0) at T = 0.3A . Here we defined A = A(T = 0). 



temperatures T > A(T) all singularities in the non-local 
conductance are smeared and essentially disappear. 

At any subgap frequency \u>\ < 2 A the zero bias non- 
local conductance decays exponentially with the distance 
between the barriers. Similarly to the case w = for 
all \oj\ < A we find IG12I ~ exp(— \x\ — x 2 1/£) as it is 
illustrated in Fig. 4. 

Non-linear response. — Finally let us turn to the 
non-linear non-local response. Choosing the bias voltage 
in the form V2 (t) = V2 + SV2 cos tot we reconstruct the 
phase 

eSV 2 . 

<P2[t) = ev2t H sinwt. 

LU 

Substituting this function into Eq. (25) and averaging 
sin[cp2(£i) — ^2^2)] over time, we obtain 



sin[cp 2 (*i) - ^2^2)] = sin[ey 2 (ii - h)} 
T 2 feSV 2 



J2 
•J I) 



eSV 2 



71=1 



OJ 



cos[nuj(ti — t 2 )] 



, (33) 



where J n {x) are the Bessel functions. Then the non-local 
differential conductance takes the form 



X 



Gu (o.r 2 + ^) +G12 (o,v 2 



uj 

nui\ 



■ (34) 



Here Gi2(0, V2) is non-local conductance defined in Eq. 
(32). We observe that the non-local differential conduc- 
tance (34) does not change its sign and for non-zero oj 
remains positive at all values of the bias voltage V2 ■ 
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Fig. 4: The absolute value of G12 at T = as a function of the 
distance between the barriers \x2 — xi\. Similarly to the case 
lu — the dependence IG12I ~ exp( — \xi — a>2|/£) is observed at 
subgap frequencies. Here we set V2 = 0, L = 20£ and x\ = 7£. 



Concluding remarks. — Our analysis of the effect 
of an external ac bias on non-local electron transport in 
multi-terminal NSN structures revealed a number of inter- 
esting features which can be tested in future experiments. 
It turns out that in the presence of an ac field exact can- 
cellation between zero-temperature EC and CAR contri- 
butions to the non-local current is lifted already in the 
lowest order in the transmissions of NS interfaces. In the 
presence of ac bias - unlike in its absence - the conduc- 
tance G12 (27) is not anymore proportional to the density 
of states in the S-metal (which vanishes at subgap ener- 
gies thus explaining the cancellation [7]). Hence, no can- 
cellation between EC and CAR contributions to G\i{uS) 
should be expected. Considering these two contributions 
separately, we observe that both EC and CAR processes 
get enhanced by external radiation, however the subtle 
balance between them is eliminated. As a result, at suffi- 
ciently small (typically subgap) frequencies ui of the exter- 
nal ac signal the non-local conductance of the NSN device 
G12 turns negative implying that CAR contribution wins 
over that of EC. Although our present analysis was car- 
ried out for diffusive structures our key conclusions should 
also apply in the (quasi-)ballistic limit. 

The behavior described above can be realized in multi- 
terminal NSN structures either under direct application of 
an ac bias or if the geometry is such that an ac Josephson 
effect can occur somewhere in the system in the presence of 
a dc voltage bias. This can be the case, e.g. if a weak link 
(such as a pinhole or constriction) is occasionally formed 
between different superconducting terminals. In this case 
ac Josephson generation may occur under dc bias effec- 
tively playing a role of an ac signal analyzed in our work. 

We also note that our results - though indirectly - could 
also be of some relevance to the effect of electron-electron 
interactions on the non-local properties of NSN devices. 
Indeed, it is well known that this effect can generally be de- 
scribed by (quantum) electromagnetic fields which medi- 
ate electron-electron interactions. In this respect photon- 
assisted violation of balance between EC and CAR (in fa- 
vor of the latter process) considered here might give a clue 



in which way Coulomb interaction could be responsible for 
the negative non-local conductance observed in recent ex- 
periments. However, a more elaborate calculation is neces- 
sary to properly account for an interplay between disorder 
and electron-electron interactions in multi-terminal NSN 
structures. This calculation will be published elsewhere. 
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